Particle simulation of vibrated gas-fluidized beds of cohesive fine powders 
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We use three-dimensional particle dynamics simulations, coupled with volume-averaged gas phase 
hydrodynamics, to study vertically vibrated gas-fluidized beds of fine, cohesive powders. The 
volume-averaged interstitial gas flow is restricted to be one-dimensional (ID). This simplified model 
captures the spontaneous development of ID traveling waves, which corresponds to bubble forma- 
tion in real fluidized beds. We use this model to probe the manner in which vibration and gas 
flow combine to influence the dynamics of cohesive particles. We find that as the gas flow rate 
increases, cyclic pressure pulsation produced by vibration becomes more and more significant than 
direct impact, and in a fully fluidized bed this pulsation is virtually the only relevant mechanism. 
We demonstrate that vibration assists fluidization by creating large tensile stresses during transient 
periods, which helps break up the cohesive assembly into agglomerates. 



o 

CO 



PACS numbers: 



I. INTRODUCTION 
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It is well known that fine, cohesive particles cannot 
be fluidized easily. Q, Q When fluidization by an up- 
flow of gas is attempted, assemblies of such particles 
tend to lift up as a plug and only form ratholes and 
cracks through which the fluid escapes. Attractive inter- 
particle forces frequently arise as a result of capillary liq- 
uid bridges or van der Waals forces. A variety of tech- 
niques to achieve smooth fluidization of such particles 
have been explored in the literature. When the cohesion 
arises from van der Waals forces, as in the case of Gel- 
dart type C particles which are typically 30 /zm or smaller 
in size, coating with hard nanoparticles [j| or an ul- 
trathin film £| can weaken the attraction between the 
bed particles, thus enabling smooth fluidization. Alter- 
nate approaches to facilitate fluidization include causing 
agglomerate break-up through a secondary supply of en- 
ergy using mechanical vibration, [ol Irj M. let Wi Ilbl ITU 
acoustic waves 13 113 or an oscillating magnetic 
field. [Ig Such approaches have been shown to be ef- 
fective even for beds of nanoparticles. ^ n the present 
study, we are concerned with some aspects of the me- 
chanics of vibrated fluidized beds of cohesive particles. 

Vibrated layers of granular materials have been stud- 
ied extensively in the literature and formation of patterns 
in shallow granular layers is now well known. |l7l ITsj 
Mixing and segregation in such vibrated beds have also 
received considerable attention. [T^. I20L l2lj The motion 
of large intruders in vibrated beds, leading to the well- 
known Brazil nut and reverse Brazil nut effects, 0, |2^ 
has also received much attention in the literature, where 
experimental measurements suggest a non-negligible in- 
fluence of the interstitial gas on the observed flow pat- 
terns. }23ll24| In general, the deeper the bed and/or the 
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smaller the bed particles, the greater the influence of the 
interstitial gas phase on the dynamics of the assembly of 
particles. 

In vibrated fluidized beds, where vibration is supple- 
mented with a fluidizing gas flow (or vice versa) , the im- 
portance of the interstitial gas is obvious as the drag 
due to the gas flow supports a substantial portion of 
the weight of the particles. Vibrated fluidized beds have 
found many applications in industrial practice (e.g., see 
an article by Squires |25|). Understanding the manner in 
which the vibration aids the fluidization process is impor- 
tant both for macroscopic analysis of vibrated fluidized 
beds and for detailed interrogation of agglomerate size 
distribution and mixing at the agglomerate and particle 
scales. 

Predicting the minimum fluidization velocity of a vi- 
brated fluidized bed (of cohesive particles) is perhaps the 
simplest quantitative, macroscopic analysis problem one 
can think of. This indeed has been a subject of many 
investigations 0, HH EH EE E3, EE HI and different 
approaches have been proposed in the literature to in- 
corporate the effect of vibration on the overall force bal- 
ance used to determine the minimum fluidization veloc- 
ity. Musters and Rietema [2(| suggested that additional 
terms need to be included in the force balance relation 
in order to account for the increased pressure drop due 
to cohesion. They included additional terms for inter- 
particle forces as well as wall friction. Liss et al. j2?J 
proposed an additional term to account for the effect 
of cohesion arising from liquid bridges. Wank et al. Q 
showed that the agglomerate size decreases as the vi- 
bration intensity increases, and studied the effect of the 
pressure on the minimum fluidization velocity. Erdesz 
and Mujumdar |2Sj developed a theory which includes 
the effect of vibration in the prediction of the minimum 
fluidization velocity; they found that the pressure drop 
decreased with increasing vibration intensity in their ex- 
periments with various particles in the range 0.15-2.75 
mm. However, Tarisin and Anuar j2?J found the oppo- 
site trend in their study of vibrofluidization of particles of 
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15 fim to 34 fiva, and others 0,0 El found no apprecia- 
ble dependence on the vibration intensity with particles 
of range 12 nm through 100 fim. 

Such conflicting reports stem from our limited under- 
standing of the manner in which the vibrating bound- 
ary interacts with the assembly of particles in the bed. 
One can readily envision at least two modes of inter- 
action: (a) When the effect of the interstitial gas on 
particle motion is of negligible importance, the vibrating 
base plate clearly imparts impulse to the particle assem- 
bly periodically only through direct collisions, (b) When 
gas-particle drag is non-negligible, a periodically varying 
pressure field can be expected to develop in the bed, and 
the base plate interacts with the bed particles indirectly 
by driving these pressure pulsations as well. The relative 
importance of these two contributions can be expected 
to depend on the bed depth, particle size, and the super- 
ficial velocity of the fluidizing gas; however, quantitative 
estimates are unavailable. 

At a more detailed level, it is readily apparent that 
vibration results in the formation of small agglomerates 
which are more amenable to fluidization, yet how such 
small agglomerates are formed is not understood. A clear 
understanding of this mechanism is an important first 
step in any effort to capture the effect of vibration on the 
dynamics of agglomerates (e.g., in a population-balance- 
type model for vibrated fluidized beds). At first glance, 
one may speculate that vibration drives vigorous colli- 
sions between agglomerates, which in turn causes their 
break-up; but it is not known if this is indeed the domi- 
nant mechanism. 

In the present study, we have examined the interac- 
tion of the vibrating base plate with a bed of particles 
through detailed simulations. We will demonstrate that 
as the fluidizing gas velocity is increased, the interaction 
with the boundary through pressure pulsation becomes 
more and more dominant and that above the minimum 
fluidization conditions, the interaction occurs almost ex- 
clusively through pressure pulsations. We will also show 
that the large tensile stress induced by the vibration is 
a more important mechanism (than vigorous collisions 
between agglomerates) in causing the break-up of the co- 
hesive particle assembly. 

The dynamic behavior of gas fluidized beds containing 
a large number of particles has been widel y e xamined 
in the literature through two-fluid models. |30l l3ll l3^ | 
This approach requires closure relations for gas-particle 
interactions and the stresses. Considerable progress has 
been made in developing and validating the closures for 
assemblies of non-cohesive particles, [3(j, LiH |32| but not 
for the case of cohesive materials. Accurate two-fluid 
model boundary conditions that capture the effect of vi- 
brating boundaries are also unavailable. Furthermore, 
the two-fluid model approach is not well suited for inves- 
tigation of the mechanism of agglomerate break-up. For 
these reasons, we model vertically vibrated gas fluidized 
beds of fine powders using a hybrid scheme, where the 
solid phase is treated as discrete spherical particles, fol- 



lowing the so-called discrete element method (DEM; or 
soft-sphere molecular dynamics) , [33| while the gas phase 
is treated as continuum, following the volume-averaged 
hydrodynamics just like in a two-fluid model. |3lL l3^ | 
This DEM-based hybrid approach was proposed by Tsuji 
et al. [34[ in their simulation of two-dimensional (2D) flu- 
idized beds, and it has been subsequently refined by oth- 
ers. [35|, Liil Recently, this model has been used to study 
fluidized beds of Geldart A particles, 37] the effect of an 
arbitrary cohesive force proportional to the particle buoy- 
ant weight, |23, and segregation in fluidized beds of 
bidisperse particles. 0] 

Even with the high speed computing available today, 
one can only simulate small systems, which are orders 
of magnitude smaller than real vibrated fluidized beds, 
through this hybrid approach. Therefore, it is essential 
that one choose suitable, idealized problems to probe the 
underlying mechanics. With this in mind, we have cho- 
sen essentially one-dimensional (ID) vibrated fluidized 
beds of particles, applying periodic boundary conditions 
in the two lateral directions. This limits the macroscopic 
dynamics to the vertical direction only, and one ob- 
tains one-dimensional traveling waves (1D-TW) instead 
of bubble-like voids observed in experiments and in fully 
three-dimensional flow simulations; nevertheless, such an 
idealized problem is, in our opinion, adequate to expose 
the manner in which the vibrating base plate interacts 
with the bed particles and how dense cohesive regions 
are broken down into smaller agglomerates. 

The rest of the paper is organized as follows. The 
DEM-based hybrid model is described in Sec. |nj and 
bubble formation, realized as 1D-TW in our geometry 
(narrow cross sectional areas), is presented in Sec. lIII Al 
The effect of cohesion and vibration on the fluidization 
is presented in Sec. IIII Bl and the manner in which the 
vibration enhances the fluidization will be discussed in 
Sec. IIII CI The pressure drop in vibrated fluidized beds 
and the mechanism through which the vibration breaks 
up cohesive assemblies will be presented in Sec. IIII Dl and 
Sec. IIII El respectively, which are followed by the conclu- 
sions in Sec. IIVI 



II. METHOD: DEM-BASED HYBRID MODEL 

Since its introduction by Cundall and Strack [3^] 
nearly three decades ago, the DEM has been success- 
fully used in modeling various particulate fl ow p roblems, 
includinghopper flows, shea ring cells, 0] rotating 
drums, |43 and oscillated layers. |44Ll45| Comprehensive 
description of this method can be found in the litera- 
ture. |1(|0,0] Below, we briefly describe the main idea 
of the DEM, and subsequently focus on how the volume- 
averaged gas phase hydrodynamics is coupled with the 
individual particle dynamics in our model. 
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A. Discrete element method 

In the DEM simulation, particles are modeled as "soft" 
sphere (the deformation is accounted for by overlaps), 
whose trajectories are computed by integrating Newton's 
equations of motion. When objects (particles or system 
boundaries) get into contact, the interaction is resolved 
by decomposing the interaction force into the normal and 
tangential directions relative to the displacement vector 
between the objects at contact, F cont = (F„,F S ), and 
the energy dissipation upon contact is characterized by 
the inelasticity and the surface friction. We use the so- 
called spring-dashpot model, with a Hookean spring. |3^| 
The objects are allowed to overlap upon contact, and the 
contact force in the normal direction F„ is determined 
by the amount of overlap A„ and the normal component 
of the relative velocity at contact v n , 



(k n A n - 7„u„)n, 



(1) 



where k n is the spring stiffness in the normal direction, j n 
is the damping coefficient, and n is the unit vector in the 
normal direction at contact, pointing from the contact 
point toward the particle center. The damping coefficient 
7„ is reiated to k n by the normal coefficient of restitution 
e (0 < e < 1), 
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where 1/m* = l/m, + 1/rrij, and i and j are indices of 
interacting particles or objects. In principle, the value 
of k n is determined by Young's modulus of the material 
under consideration. However, unless stated otherwise, 
we use a much smaller value for k n , compared to the one 
computed based on the usual range of Young's modu- 
lus. If the main results are not qualitatively different, 
it is favorable to use a smaller value of k n , because the 

— 1/2 

collision duration time in DEM scales with k n , which 
determines the integration time-step size required to ac- 
curately resolve the interaction during the contact. This 
is a well-known issue in DEM simulations. 0, ^3 We 
varied k n over three orders of magnitude and verified 
that the main results do not depend sensitively on the 
choice of k n , even though actual contact force between 
the objects certainly depend on the value of k n . For in- 
stance, the results for two k n values differing by a factor 
of 10 will be presented in Figs. 003 and El 

The interaction in the tangential direction is modeled 
by a "spring and slider", and the contact force is given 
by: 



F s = -sign(v s ) x min(fc t A s ,^|F„|)s, 



(3) 



where v s = v s ■ s is the tangential component of the rela- 
tive velocity at contact; v s = n x (vy x n); s is the unit 
vector in the tangent plane collinear with the component 
of the relative velocity at contact; k t is the tangential 
spring stiffness that is related to k n by the Poisson's ra- 
tio of the material vp [kt — 2k n (l — vp)j{2 — vp)}\ and A s 



is the magnitude of tangential displacement from the ini- 
tial contact. The magnitude of the total tangential force 
is limited by the Coulomb frictional force /i|F„|, where /i 
is the coefficient of friction. More sophisticated and/or 
realistic interaction models, such as that of Walton and 
Braun's 0] or a Hertzian spring-dashpot model, jH3 
may also be used. However, we choose a simple Hookean 
spring-dashpot and spring-slider model, as it has been 
shown to successfully reproduce many experimental ob- 
servations, 0, 0j C3 an d it is computationally more 
tractable than others. 

Among different inter-particle forces, other than due 
to contact, we consider only cohesion arising from van 
der Waals force. In principle, the cohesion can depend 
on the particle characteristics, such as polarizability, par- 
ticle size, and asperity. However, we adopt a simple 
formula by Hamaker, as we aim to bring out the effect 
of cohesion on the fluidization behavior, rather than to 
validate different cohesion models. Particulate flows in 
industry often consist of particles with a wide range of 
sizes and shapes; however, we seek a better understand- 
ing of simple systems consisting of monodisperse spheres, 
which are well characterized by a small set of parameters. 
The cohesive van der Waals force between two spheres of 
radii and rj can be expressed as, |5ll | 



A 2r l r j {s + r i +r 3 ) 
' 3 [s(s + 2 n + 2r )} 2 

s{s + 2r l + 2r j ) 



x 



(s + n + rj) 2 - (n - rj) 2 

A r . 

-n (tor r = r% = rj 



n 



12 s 2 



- 1 

r-i = rj and s<r), (4) 



where A is the Hamaker constant, and s is the inter- 
surface distance. As the original formula is a rapidly 
decreasing function of s, further simplification using the 
assumption of s < r has been made. This model has 
a singularity at contact. In order to avoid this artifact, 
we introduce a widely accepted minimum cut-off value 
for the inter-surface distance of 0.4 nm (= 5*), which 
corresponds to the inter-molecular center-to-center dis- 
tance, [j^ In what follows, the level of cohesion is repre- 
sented by the cohesive Bond number Bo, which is defined 
as the ratio of the maximum cohesive force (at the mini- 
mum cut-off separation 6*) to the particle weight. Other 
types of cohesion can be readily accounted for in DEM- 
based models. IMElls! 



B. Coupling with gas phase hydrodynamics 

The dynamics of individual particles is coupled with 
the volume-averaged gas phase hydrodynamics. |34| In 
this hybrid model, the equations of motion for individual 
particles have two additional terms (compared to tradi- 
tional DEM modeling particles in vacuum) arising from 
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the presence of the gas phase: 



' dt 



m p g e ff+F c 



V. 



m<t>) (u fl 



where m p and v p 



P )-V p Vp, 
(5) 

are individual particle mass and ve- 
locity, respectively. The first term on the right hand 
side represents the body force due to gravity; g e ff is 
the effective gravitational acceleration in the reference 
frame where equations are integrated. For non-vibrated 
beds, g e ff is simply the gravitational acceleration g. 
When the bed is subject to a single frequency oscilla- 
tion, the equations are integrated in the vibrated frame, 
and g eff = g[l + rsin(2ir/t)], where T = A p (2nf) 2 /g 
is the maximum acceleration of the base plate (distrib- 
utor) non-dimensionalized by the gravitational accelera- 
tion g = |g|, A p is the oscillation amplitude, and / is 
the oscillation frequency. We assume the oscillating base 
plate is made of the same materials as the particles (the 
same values for e and //), and that the mass of the plate 
is infinitely large compared to that of an individual par- 
ticle. The second term and the third term represent the 
aforementioned contact force and van der Waals force, 
respectively. The total force acting on the particles due 
to the fluid is commonly partitioned into the local drag 
part and the effective buoyant part, as was done here 
(see e.g., an article by Ye et al. :37j|) : The fourth term 
accounts for the drag force, and the last term accounts 
for the contribution of the gradually varying part of the 
pressure field, where V p is the volume of each particle; <\> 
and u 9 are volume-averaged solid phase volume fraction 
and gas phase velocity, respectively; (3 is the inter-phase 
momentum transfer coefficient [3l| ; and p is the gas phase 
pressure. 

In general, the gas phase quantities are obtained by 
simultaneously integrating the coarse-grained mass and 
momentum balance equations. We assume the gas phase 
to be incompressible, which will be validated later (see 
Sec, imrjjl . The addition of continuity equations for the 
gas phase and solid phase reads 



V-[(l 



0, 



(6) 



and a reduced momentum balance equation for the gas 
phase, based on generalized Darcy's law, is given as 



= -(!-<£) Vp + /3(<£) (u s -u 3 ), 



(7) 



where u s is the coarse-grained solid phase velocity. 
Coarse-grained variables are considered only on grids 
where the continuum balance equations are solved. Note 
that the solid phase continuum (or coarse-grained) vari- 
ables are explicitly available in the course of DEM com- 
putation. 



C. Beds of narrow cross sectional area 

Only beds of narrow cross sectional areas will be con- 
sidered, and the volume-averaged gas phase (hence the 



solid phase coarse-grained variables as well) is assumed to 
be ID. However, the solid phase is maintained to be 3D, 
as the way particles pack and collide in lower dimensions 
are considerably different from those in realistic 3D cases. 
Our assumption allows us to consider relatively deep beds 
(through inexpensive computational effort) and to bring 
out the basic physics of more complicated dynamics in 
higher dimensions. 

Solid phase coarse-grained variables at ID discrete grid 
points are computed by distributing the particle mass 
and momenta to the nearest two grid points using a halo 
function h that continuously decreases to zero around the 
particle; 



h(z - z ) 



1 — \z — zq\/ Az 




for \z — zq\ < Az, 
otherwise, 



. (8) 

where z is the particle position in the vertical direction, 
zq is that of a neighboring grid point, and Az is the grid 
spacing. It is readily seen that h has the property that the 
particle quantities are distributed to the two nearby grid 
points, inversely proportional to the distance to the grid 
point. The coarse-grained variables, the number density 
n and u s , on the grids are then defined simply as 

N 



n ( z a) = ^2h(zi - z Q ), (9) 
»=1 
N 

i(z )u s (z ) = h(zj - Zo)vp,j. (10) 



where zi and zq are the ith particle location and nearby 
grid location, respectively. 

The assumption of the gas phase to be ID facili- 
tates further mathematical simplifications of the above 
particle-gas interaction formulation. In ID continuum 
cases, Eq. © can be integrated 



(1 - (j))U g + (j)U S = U S 



(11) 



where XJ S is the superficial gas flow velocity. Strictly 
speaking, in a vibrated fluidized bed, U s may also vary 
periodically. The extent of its variation will depend on 
the dynamics of the gas in the plenum and the flow re- 
sistance offered by distributor (base plate). One can 
show that the temporal variation of U s will be small 
for a highly resistive distributor plate, which we as- 
sume. Thus, in modeling of both non-vibrated and vi- 
brated fluidized beds, U s will be considered to be a time- 
independent parameter. 

After some manipulation, Eq. (JSJ) can be rewritten as 
follows: 



~dt 



V, 



Ipgeff 



K - V p) - (1 _ 6)2 ~ U * 



(12) 



Note that the presence of the gas phase is realized as 
additional terms involving coarse-grained variables, in- 
stead of separate continuum equations to be integrated 
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simultaneously. In the course of integration, <f> and u s 
in Eq. I112H need to be evaluated at the particle location, 
not at the grid points. We evaluate them by linearly 
interpolating those values at the neighboring grid points. 

For the inter-phase momentum transfer coefficient /3, 
we use an expression proposed by Wen and Yu [55j : 







3 - PgH 1 



</>)l u s 



-(1 



-2.65 



(13) 



where Co is the drag coefficient, p g is the gas phase mass 
density, and d p is the particle diameter. The drag coef- 
ficient proposed by Rowe [5(j is employed in our model: 



U =0.041 U =0.082 U =0.12 
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time time time 



C D = 
where 



#- (l + 0.15i?e' 



24 
Re 

0.44 



0.687\ 



Re g < 1000, 
Re g > 1000, 



Re a 



(1 - (f))pgd p \u.g - U S 

Hg 



(14) 



(15) 



and p g is the gas phase viscosity. As we consider fine 
powders, Re g is generally small, and we use the assump- 
tion Re g <C 1, which further simplifies B: 



m = i8^(i 



-2.65 



(16) 



Casting Eqs. Q12II and (|16[1 in a dimensionless form, 
using p s , dp, y/gdp, y/d p /g as characteristic density, 
length, velocity, and time, one obtains the following non- 
dimensional groups (arrows indicate changes in the nota- 
tion from dimensional variables to non-dimensional vari- 
ables that will be used henceforth): 



Psgdl ' 



6 = T> 
dp 



Bo 
St 



A 



4irp s gd p S* 2 
Psy Up 

Hg 



spring stiffness 
superficial gas flow rate 
scaled minimum separation distance 
cohesive Bond number 

Stokes number 



together with non-dimensional parameters, namely T, 
f *~ fy/dp79, e, /it, and v P . 



III. RESULTS AND DISCUSSION 

We simulate both non-vibrated and vibrated gas- 
fluidized beds of non-cohesive or cohesive particles (0 < 
Bo < 50). We consider beds of narrow square-shaped 
cross sectional area 5d p x 5d p or 10(ip x 10d p with static 
depth Hq ~ 100dp and Hq ~ 200d p (which consist of 



FIG. 1: Spatiotemporal plots of conventional fluidized beds 
of non-cohesive particles (Bo = 0.0, St = 55) at three differ- 
ent superficial gas flow rates. Gray scale represents the local 
volume fraction of particles; both the regions of the minimum 
value ((f) ~ 0.33) in the gray scale bar and completely void 
regions ((f> = 0) above the bed top surface are shown in the 
same white color. 



3000 and 6000 particles respectively, in beds of 5d p x 5d p 
cross section) . Periodic boundary conditions are imposed 
in both lateral directions, in order to avoid strong side 
wall effects in beds of such a small aspect ratio. We check 
that our results do not sensitively depend on a particular 
choice of cross sectional area or the depth of the bed (e.g., 
see Fig- EJ ■ We mostly use a bed of ~ 5d p x 5d p x 100d p 
in the following computations, unless otherwise stated. 
We used d p , 1.5d p , and 2d p for the grid spacing Az. The 
detailed profiles of the coarse-grained variables slightly 
depend on the choice of Az (the bigger the grid size is, 
the smoother the variables are, as one can readily ex- 
pect), but the main results remain the same, unless the 
grid size is too large; we set Az = 1.5d p in all the results 
presented here. In the following, all the quantities will 
be shown in non-dimensional form. 



A. Bubbling and one-dimensional traveling waves 

We start by considering conventional (non-vibrated) 
gas fluidized beds of non-cohesive particles to make cer- 
tain that our model captures basic experimental obser- 
vations. We first estimate the minimum fluidization ve- 
locity through a simulation of quasi-static increase in the 
gas flow rate (from zero) and the measurement of the 
pressure drop across the bed, which yields U m f w 0.022. 
This estimate is slightly smaller than what we can com- 
pute using the force balance relation and the approximate 
formula of Wen and Yu, U m f « 0.023 [57|; in this calcu- 
lation, we used (f> m f = 0.652, which is measured from the 
bulk of the bed. When the gas flow rate exceeds U m f, 
the bed in our model starts to expand inhomogeneously, 
and forms alternating bands of plugs and voids. Particles 
located at the bottom of one plug "rain down" through a 
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FIG. 2: Traveling wave speed (filled circles) in non-vibrated 
gas-fluidized beds (F = 0), for various levels of cohesion. U s ~ 
0.12. Beds of larger Bo's (>~ 6) do not get fluidized; rather 
the whole bed moves up as a plug. Insets are spatiotemporal 
plots of moderately cohesive bed that is fluidized (Bo — 2) 
and more cohesive bed (Bo = 8) that is not fluidized and 
rises as a plug at this flow rate. Triangles: The wave speeds 
obtained in beds of ~ 5d p x 5d p x 200d p . Crosses: The wave 
speeds obtained in beds of • 



10d p x Wd p x lOOdp. 



void region and accumulate at the top of the lower plug, 
causing the void regions to rise to the top in a periodic 
fashion (Fig. This phenomenon corresponds to the 
formation of a periodic train of bubble-like voids in real 
fluidized beds, which appears as 1D-TW in the narrow 
beds we consider. These waves represent the first stage 
in the bifurcation hierarchy leading to various inhomo- 
geneous structures in hi ghe r-dimensional (i.e. 2D or 3D) 
fluidized beds. [58L Ifj9l l60j] The secondary bifurcations 
which occur in real fluidized beds are suppressed in our 
simulations, which are forced to retain the ID character. 
As the gas flow rate increases, both the wave speed and 
amplitude increase (Fig.^l. In the subsequent Sections, 
we will use the 1D-TW as an indicator that characterizes 
the fluidizability of the bed. 



B. Cohesion, vibration, and fluidization 

When the gas flow rate is well above U m f, the bed 
exhibits clearly visible 1D-TW [for instance, see the cases 
of U s = 0.082 and 0.12 in Fig. [T]. In this Section, we 
examine the influence of the cohesion (between particles) 
on the fluidizability of a bed, and explore how mechanical 
vibration facilitates the fluidization of beds of cohesive 
particles. 

Figure [21 shows the effect of Bo on the wave speed, 
where we have kept St ~ 55 and U s — 0.12 (the same val- 
ues as in the last panel of Fig.^J. As Bo increases, the 
wave gradually slows down, and eventually disappears. 
At Bo = 6, the wave travels intermittently, remaining 
stationary for some time and traveling at other times, 
and the wave speed during the non-stationary phase is 



FIG. 3: Traveling wave speed in beds of highly cohesive 
powders (Bo — 20) subject to vibrations with various levels 
of r. The traveling wave reappears for T > 2. Insets are 
spatiotemporal plots of three different cases (r = 3.0, 6.0, 10), 
shown in the vibrating plate frame. / = 0.16; U s = 0.12. 



shown in Fig. [21 For Bo > 7, the whole bed rises up 
as a plug at this flow rate, which is consistent with the 
well-known experimental observations in narrow beds of 
strongly cohesive particles (see a review article by Sun- 
daresan [6(| and references therein). Further increase of 
U s only slightly improves the fluidizability of the bed, 
confirming that beds of highly cohesive particles cannot 
be fluidized by simply increasing the gas velocity. 

The results shown in Fig. [21 are obtained with the as- 
sumption that the interaction between the particle and 
the base plate is non-cohesive. We check that the ad- 
hesion at the base plate does not make any difference in 
the above results, as well as all the main results presented 
here. The only noticeable difference is that strongly co- 
hesive particles at the bottom of a bed get stuck to the 
plate for some time (during a cycle; in vibrated beds), 
depending on the oscillation parameters and the superfi- 
cial gas flow velocity. Detailed comparison between the 
two cases (with and without adhesion) is shown later in 
Fig. El 

Now we subject beds of even more cohesive particles 
(Bo = 20) to mechanical vibration of a single frequency 
sinusoidal oscillation in the direction of gravity. When 
the vibration intensity is strong enough (when T exceeds 
a certain value), even these highly cohesive beds get flu- 
idized in the sense that 1D-TW reappears (Fig. 0). At a 
fixed gas flow rate (we assume it to be time- independent; 
see the discussion in Sec. Ill C(l . the wavelength apparently 
increases with T; however, the wave speed remains nearly 
the same (Fig. |3J). We define the critical Bond number 
Bo c as the maximum value of Bo for which the bed is 
fluidizable (generating 1D-TW) at given set of oscilla- 
tion parameters, and compute it as functions of T and /, 
using bisection-type search along the .Bo-axis. We find 
that Bo c increases almost linearly with T (Fig. , but it 
only weakly depends on / (Fig.EJ. The frequency range 
shown in this figure corresponds to the usual operation 
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scaled plate trajectory 



FIG. 4: The critical values of Bo, above which the bed does 
not get fluidized at given vibration parameters, as a function 
of r. Three different levels of Us are used for F = 3.0 and 
6.0. / = 0.16. A solid line, the least square fit for the cases 
of Us = 0.12, is included to guide the eye. 



range between 10 Hz and 100 Hz, for the particle size of 
d p = 50 /im. In this range of frequencies, Bo c is virtu- 
ally independent of the frequency. As the gas flow rate 
increases, Bo c slightly increases at fixed values of T and 
/ (see the cases for T = 3 and 6 in Fig.BJ. 



C. The role of vibration: Direct impact vs. 
pressure pulsation 

In this Section, we discuss the manner in which the 
vibrating base plate or the distributor, interacts with the 
bed material. 

In the absence of gas, the kinetic energy of individual 
particles in vibrated beds is obtained only from direct im- 
pact with the plate, and dissipated through interaction 
between particles. In shallow beds (Hg <^ 20d p ), the 
fluctuating kinetic energy (granular temperature) dissi- 
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FIG. 5: The critical values of Bo as a function of /, for two 
different levels of F. U s = 0.12. 



FIG. 6: The bottom surface and the average vertical velocity 
of vibrated beds (averaged over all the particles in the bed) 
during a cycle, shown in the vibrating frame, obtained for 
three different gas flow velocities. Bo = 20; F = 3.0; / = 0.16. 
A thin solid line, an arbitrarily rescaled plate trajectory in the 
laboratory frame, is drawn to represent the phase angle during 
a cycle. 



pates so quickly through collisions that the whole bed can 
be well approximated by one solid body. The dynamics of 
the center of mass of such a bed can be described by that 
of single perfectly inelastic ball on a vibrating plate. [(3l| 
Vibrated shallow beds in vacuum undergo period dou- 
bling bifurcations as T is varied. [63] Such layers of a 
large aspect ratio form various spatiotemporal standing 
wave patterns. 0, 0] However, temporal dynamics of 
vibrated deep beds (of about 20 particle deep or more) 
in vacuum are not commensurate with the oscillation fre- 
quency, and exhibit more complicated non-periodic be- 
havior. [r33| 

When the gas phase effects are accounted for, the gas 
drag causes the dynamics of vibrated deep beds to de- 
viate from those in vacuum, and the deviation is more 
pronounced for smaller particles (because the gas drag is 
larger). For deep beds of fine powders considered here, 
the presence of gas phase modifies the bed dynamics so 
that the temporal dynamics become periodic. In the ab- 
sence of a net flow (U s = 0) the bed lifts off from the 
plate only slightly (even smaller than the particle size) 
during a fraction of a cycle, and the bed impacts the plate 
later in the same cycle (dot-dashed lines in Fig.[f)|). This 
periodic behavior has the same periodicity as the plate 
oscillation. Even for higher values of T, period doubling 
phenomenon, which occurs in vibrated shallow layers in 
vacuum, is not observed. Note that, in the early phase 
of the oscillation cycle shown in Fig. 0(0 < t/T < 0.25), 
the plate moves downwards (and so does the base plate), 
and yet the bottom surface of the bed is approaching the 
distributor plate, i.e. the bed is descending faster than 
the plate. The bed hits the plate, stays in contact for a 
duration of time, and then detaches from it. 

When the gas flow is turned on and its rate gets in- 
creased, the velocity of the bed (relative to the base plate) 
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FIG. 7: (Top) Normal stress tensor on the base plate, scaled 
by the bed weight per unit cross sectional area (Emj), is 
shown during a cycle for different superficial gas flow rates. 
Bo = 20; T = 3.0; / = 0.16; k„ = 2.0 x 10 5 . Solid lines 
are obtained by neglecting adhesion at the base plate (Bo p , 
the Bo at the plate, is 0). Dashed lines are obtained when 
the adhesion is accounted for (Bo p = 40), for two cases 
(Us = 0.0082, 0.024), which are nearly the same as the non- 
adhesive cases. A dot-dashed line, the plate trajectory in an 
arbitrary unit, is drawn to represent the phase angle during a 
cycle. (Bottom) The scaled pressure drop for the same cases 
as in the top panel. Inset shows the sum of the pressure 
drop and the normal stress at the base plate for three cases 
(Us = 0.0,0.016,0.033), all of which virtually coincide with 
the reduced effective gravity (g e ff/g) in the vibrating frame 
[sinusoidal curve; 1 + rsin(27r/t)]. Small excess amounts, ap- 
pearing as small peaks on top of a sinusoidal curve, arise from 
the transient restoring force of the soft particles at the impact, 
which gradually disappears as the plate impacts the bed more 
gently with increasing U s (there are no apparent additional 
peaks for U a = 0.033). 



during its short flight increases, and the bed lifts off fur- 
ther from the plate (Fig.EJ). The upward gas flow resists 
downward motion of the bed, hence not only the dura- 
tion of direct impact but also its magnitude (strength) 
gradually decreases (top panel in Fig. |JJ , as the gas flow 
rate increases. When the adhesion at the base plate is 
accounted for (dashed lines in Fig. [7J, the plate experi- 
ences some force even when the bulk of the bed is in flight 
(t/T <~ 0.2 or >^ 0.6), as there are a small number of 
particles stuck to the plate. Other than this, compared to 
the case when the adhesion is neglected (solid lines), no 
difference is observed. The bed eventually hardly touches 



the plate at some gas flow rate, above which direct impact 
remains minimal. It will be shown in the next Section 
that the minimum gas flow rate at which the direct im- 
pact virtually vanishes is for all practical purposes the 
same as the minimum fluidization velocity in vibrated 
beds. 



D. Pressure drop in vibrated beds 

The pressure drop across vibrated beds oscillates with 
the same periodicity as the plate oscillation (bottom 
panel in Fig. 01. The pressure drop increases (decreases) 
when the base plate moves down (up); see Fig. [5] and the 
bottom panel in Fig. El for t/T <^ 0.5. At first glance, 
this seems to be counter-intuitive, as the pressure drop 
increases when the plate is "moving away" . However, it 
should be noted that the change in the pressure drop is 
determined by the change in the gap between the bed and 
the base plate (i.e. the relative motion with respect to 
the plate), not by the absolute motion of the base plate 
in the laboratory frame. As noted in the previous Sec- 
tion, the bed approaches the plate during the phase of 
the oscillation cycle when the plate is moving down from 
its mean position. 

As soon as impact occurs, the pressure drop begins to 
decrease rapidly, even below zero. While the bed is mov- 
ing away from the plate after the take-off (see t/T >^ 0.6 
in Fig. 0), the pressure drop continues to decrease. Dur- 
ing this time, a region of lower pressure is being created 
in the gap between the bed and the plate (bottom panel 
in Fig. [3). Our results are generally consistent with the 
experimental measurements by Thomas et dl. |64| : how- 
ever, direct comparison with their data is not possible 
because of differences in systems and particle properties. 
As the gas flow rate increases, the abrupt drop at the im- 
pact gets smaller until it virtually vanishes (and so does 
a sudden increase in the stress at the plate a p i a te), and 
the pressure drop curve during a cycle approaches nearly 
the same, asymptotic sinusoidal curve. 

When the vertical component of the stress at the plate 
is combined with the pressure drop for each case, the 
resulting curves for different gas flow rates virtually co- 
incide with the same sinusoidal curve (inset in bottom 
panel of Fig. 0. This can be understood from the fol- 
lowing force balance relation in the direction of gravity: 

J2 m 9 J2 m 9 9 

where ^2 mg is the weight of the bed per unit cross sec- 
tional area, and the restoring force due to the softness 
of the particles (i.e., the spring stiffness) is neglected. 
The above relation holds at every moment during a cy- 
cle, and the asymptotic common sinusoidal curve in the 
inset is identified to be 1 + T sm(2n ft) . It corresponds to 
the effective gravity in the vibrating plate frame g e ff/g, 
the right hand side of Eq. 1)17(1. Using the fact that the 
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FIG. 9: The average pressure drop during a cycle, scaled 
by the bed weight per unit cross sectional area, obtained 
for different cases of vibrated fluidized beds, as a function 
of Us- The pressure drop in conventional fluidized beds of 
non-cohesive particles (diamonds) is included for comparison. 
Filled circles are for the case of Bo = 20, F = 3.0, but 10 times 
larger value of spring stiffness (k n = 2.0 x 10 6 ), compared to 
all the other cases (k„ = 2.0 x 10 5 ), was used. 



FIG. 8: The same quantities measured for the same cases as 
in Fig. |7| except that the spring stiffness k n and k 3 have been 
increased by a factor of 10 (e.g., k„ — 2.0 x 10 ). 

pressure drop across the bed is limited by (1 + T) ^2 mg 
during a cycle, we can examine the validity of the in- 
compressibility assumption for the gas phase that we 
use: By making use of the equation of state for an ideal 
gas, one can show that the ratio of the change in gas 
phase density Ap g to its reference value p re f satisfies 
Ap g /p ref ~ Ap/P atm < (l + T)Y,mg/P a tm, where P atm 
is the atmospheric pressure. For the beds of fine powders 
considered in our study, J2mg/P atm ~ O(10~ 5 ), hence 
Ap/ p re f -c 1 and the assumption of the incompressibility 
is valid at every moment during a cycle. 

In order to test the sensitivity of the results to the 
spring stiffness, we repeated the calculations shown in 
Fig. with an order of magnitude larger value of k n . 
The results are presented in Fig. El The detailed behavior 
during a cycle surely depends on the value of the spring 
stiffness; the beds of softer particles in Fig. [7| noticeably 
further compress and expand during the impact, and the 
pressure drop keeps decreasing even below zero during 
the impact. This effect diminishes when a larger value 
is used for k n ; in spite of the quantitative changes that 
are readily apparent, the results shown in Figs. [7| and El 
are qualitatively similar. Importantly, in a fully fluidized 
state, pressure pulsation is the only relevant mechanism 
in both cases, and the value of k n becomes irrelevant. 

Viewing the pressure pulsation at the bottom plate 
as a forcing set up by the plate, one can inquire about 
the speed at which this pulsation propagates upwards 
through the bed. As we have taken the gas phase to 



be incompressible, the pulse propagates nearly instan- 
taneously. Thus, at each time instant, the gas pres- 
sure decreases essentially monotonically as one moves up 
through the bed (except for the small periodic variation 
associated with the voidage waves). If one allows the 
gas to be compressible, then the pulsation travels at the 
speed of sound; sound speed through gas-fluidized beds 
is considerably smaller than that through a column of 
gas. [63 If the time required for the propagation of the 
pulsation through the bed is commensurate with the pe- 
riod of the plate oscillation, resonance can set in; how- 
ever, such resonance is suppressed in the present study as 
we have treated the gas as incompressible. In any case, 
in the relatively shallow beds that we consider, resonance 
is not expected to be a significant effect. 

The pressure drop averaged throughout a cycle as a 
function of gas flow rate manifests a linear increase up 
to the constant plateau region (Fig. |3J), which is quali- 
tatively the same as in conventional fluidized beds. The 
minimum velocity when the direct impact virtually does 
not occur is essentially the same as the minimum fluidiza- 
tion velocity in vibrated beds; only the pressure pulsation 
is a relevant mechanism for a fluidized state of vibrated 
beds. Note that the minimum fluidization velocities for 
the vibrated beds of cohesive particles are larger than 
that for a non-vibrated bed of non-cohesive particles of 
the same size. This can be interpreted as the increment in 
effective particle size, which is understandable, because 
the fluidized entities in vibrated beds of cohesive particles 
are agglomerates, not individual particles. The average 
volume fraction at the minimum fluidization is smaller 
(e.g. <p m j — 0.631 for the case of Bo = 40), compared 
to what is observed in a bed of non-cohesive particles 
(0.652); cohesive beds tend to pack more loosely. In the 
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FIG. 10: [Color online: Particles are color-coded according 
to the vertical velocity, increasing from blue (moving down- 
wards) to green to red (moving upwards)] Snapshots of beds of 
(a) non-cohesive particles (Bo — 0) and (b) cohesive particles 
(Bo = 4) , obtained from particles in a box of height L = 18d p 



and cross sectional area 10d„ 



10d B 



= 0.40). Periodic 



boundary conditions are imposed in all three directions (see 
text). Particles rain down from the upper plug to the lower 
plug, while void regions in the middle rise up, which form trav- 
eling voidage waves. Note that cohesive particles form strings 
or agglomerates while they rain down. Full animations are 
available from http://multiphase.princeton.edu/ICE_2005 



fully fiuidized state, the pressure drop exhibits a plateau, 
which approximately equals the weight of the bed per 
unit cross sectional area (Fig. El, as in non- vibrated beds. 
There is no clear consensus on this issue in experimen- 
tal studies. Tasirin and Anuar 01 found that the pres- 
sure drop increases as the vibration intensity T increases, 
Erdesz and Mujumdar ^ observed the opposite trend, 
and Marring et al. Mawatari et al. .9j, and Nam et 
al. |ll| observed constant plateau pressure drop equaled 
with the bed weight per unit cross sectional area in high 
gas flow rates. Only the latter is consistent with our re- 
sults, which can be explained by the simple force balance 
argument in Eq. I|17ll . accounting for the fact that the 
direct impact is negligibly small in a fiuidized state. 



E. Break-up of cohesive assembly by pressure 
pulsation 

We seek to understand how vibration facilitates the 
break-up of cohesive assemblies into agglomerates and 
maintains the propagation of the wave in a fiuidized state. 
We analyze the profile of continuum variables, including 
the granular temperature T and the solid phase stress 
tensor a (or the solid phase pressure), across the travel- 
ing wave during a cycle. As the waves in a bed of finite 
depth that we have considered thus far are not perfectly 
periodic, we consider an alternate, idealized geometry, 
where a wave is fully developed in a small periodic box (in 
all of three directions) of height L that is commensurate 
with the wavelength obtained in the vibrated fiuidized 



bed simulations described above. Note that the direct 
impact does not play an important role for a fully fiu- 
idized bed fSec. HIlCl only the pressure pulsation does), 
and the vibrating plate does not have to be considered 
in such a case. For a fully fiuidized state, the weight of 
the bed per unit cross sectional area is supported by the 
pressure drop: 



p|_=o -Pi 



Pp9eff(f>avgL, 



(18) 



where 



'avg 



is the average volume fraction. 



Comparison between fully fiuidized states of cohe- 
sive beds and non-cohesive beds on microscopic level 
reveals that cohesive particles form strings of particles 
or agglomerates while they rain down through void re- 
gions, whereas non-cohesive particles come down indi- 
vidually (Fig. 1101 full animations are available from 
http://multiphase.princeton.edu/ICE_2005). As one can 
readily see, this effect arises from the attractive force 
between cohesive particles, which can be well character- 
ized by tensile stress on a continuum level. We compute 
the stress tensor, consisting of a kinetic or dynamic part 
and a virial or static part, using the following microscopic 
relation I6rl: 



1 

V 



5> 

c€V 



(19) 



where = v,— < V; > is the fluctuating velocity of 
the ith particle, ® is the dyadic tensor product, f c is the 
interacting force between contacting particles 1 and 2, 
and l c = ri — Y2 is the displacement vector between the 
centers of particles under consideration. The second term 
is summed over all the contacts in the averaging volume 
V. 

Now we consider the continuum level interpretation of 
vibrated fiuidized beds of highly cohesive particles in a 
fully periodic box during an oscillation cycle (Fig. fTT|> . 
During a cycle, the wave oscillates up and down with a 
net upward motion, as can be seen in cases of Fig- EI We 
compute the solid phase continuum variables including 
the pressure, the trace of the stress tensor per dimension 
Tt(a)/D. We compute them in the co-traveling frame 
(with the wave's net motion), averaged over 100 cycles 
for the purpose of variance reduction. While the pressure 
pulsation cyclically varies during a cycle, so do both the 
granular temperature and the stress tensor. 

Fieure lTTl shows the variation of solid phase continuum 
variables at six different times during an oscillation cycle. 
The inset in each panel shows the position of a hypothet- 
ical, oscillating base plate at the instant the profiles of 
the continuum variables are shown. The plate trajectory 
sketched in the insets can be compared to that shown in 
the top panel of Fig. \7\ The corresponding pressure drop 
across the periodic domain can be found from the cycli- 
cally varying pressure drop profile in the bottom panel of 
Fig. [71 which is nearly out of phase with the plate posi- 
tion. In Fig. ^2 ( a ) > 9ef f/g and the scaled pressure drop 
are nearly unity; they are larger in Fig. 1111 (b); attain 
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FIG. 11: Solid phase continuum variables of a highly cohesive bed (Bo = 20, 4>avg = 0.40) during a cycle in the vibrating frame 
(r = 3.0), obtained in the same geometry as in Fig. 1101 Insets represent the phase angle during a cycle in terms of the vertical 
location of the (virtually existing) base plate (indicated by dots). The void region in the middle oscillates up and down, and 
has an upward net motion [indicated by a gray arrow in panel (a)]. The solid phase stress during a cycle can be summarized 
by the following stages: (a) The pressure pulsation is weak and the particle phase pressure is negligibly small throughout the 
bed. (b) As the pressure pulsation increases, so do the magnitude of the normal stress both in compressive (in lower plug) 
and tensile regions (in upper plug) and granular temperature, (c) until they reach their maximum values; around this time the 
tensile stress reaches the tensile strength of the material, which breaks up the cohesive assembly into agglomerates, (d) Then 
the stress relaxes, (e) and (f) Relatively small compressive stress builds up at the bottom of the upper plug, but this is not 
relevant for the break-up of the assembly. The location of the peak tensile stress is indicated by a broken arrow in panel (c). 
A full animation during two cycles is available from http://multiphase.princeton.edu/ICE_2005 



their largest values near the instant shown in Fig. 1111(c): 
relax back towards unity in Fig. 1111 (d): and then they de- 
crease in Fig. 1111 (c) until they reach the minimum values 
in Fig. El (f)- 

When the scaled pressure drop is nearly unity [Fig. 1111 
(a)], the solid phase pressure is negligibly small through- 
out the bed. As the pressure drop increases until it 
reaches its maximum value [Figs. ^2 (b) and (c)], the 
granular temperature and compressive stress in the lower 
plug and tensile stress in the upper plug also increase sig- 
nificantly until they reach their maximum magnitudes. 
The large granular temperature occurs at the lower plug, 
where the particles or agglomerates get accumulated. It 
arises from vigorous collisions among raining down par- 
ticles or agglomerates and the lower plug. Note that the 
volume fraction in this region is still low (<f> <~ 0.3). 
A rough estimate for the magnitude of the stress using 
the values at the maximum temperature region (~ nT = 
6/(irdp)<fiT ~ 0.2), where the most vigorous collisions 
occur, shows that it is still far smaller than the tensile 
strength of the cohesive assembly (~ 1) that is estimated 
below; here n is the number density. In order for these 
vigorous collisions to contribute to the break-up of the 
cohesive assembly, the magnitude of the stress formed by 
the collisions has to be comparable or larger than the 
tensile strength of the material. However, the stresses 



formed by the collisions are not strong enough to break 
up the assembly, and irrelevant for the assembly break- 
up. Rather, it is the increased tensile stress in the upper 
plug that breaks up the assembly into agglomerates and 
maintains the wave propagation. Particles in the upper 
plug can split off from the assembly and rain down, be- 
cause the increased tensile stress becomes large enough 
to reach the tensile strength of the cohesive assembly. It 
occurs at a time around that shown in Fig. ^2 ( c ) ■ 

As the magnitude of the tensile stress cannot exceed 
the strength of the assembly, we estimate the strength 
of the fluidized bed by measuring the maximum tensile 
stress. We compare such obtained tensile strength of the 
fluidized bed with the prediction of Rumpf 's model, [§3 



7T d p 



(20) 



where at is the tensile strength, e is the porosity (= 1 — <j>), 
k is the coordination number, and Ft is the cohesive con- 
tact force. We find that our measurement (a t ~ 0.8) 
is about an order of magnitude smaller than the predic- 
tion (at ~ 6). This discrepancy is understandable, be- 
cause the cohesive assembly in a fluidized bed breaks up 
through the weakest linkage, as opposed to all directions 
as is assumed in the Rumpf 's model. 
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Later in the cycle, when pressure drop decreases and 
approaches the minimum value [Figs. ^2 (d) through (f)], 
the particle phase pressure is generally small in the bed. 
When the pressure drop is negative, relatively small com- 
pressive stress builds up at the bottom of the upper plug 
as the gas in the void region pushes up; this is obviously 
irrelevant for the break-up. 



of realistic industrial scales, using the current approach, 
is not yet feasible. 
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IV. SUMMARY 

We have used a particle dynamics-based model for vi- 
brated gas-fluidized beds of fine powders to study how 
the vibration facilitates the fluidization of beds of cohe- 
sive powders. We have demonstrated that, as the gas 
flow rate increases, the direct impact from the plate de- 
creases and the pressure pulsation becomes more domi- 
nant. In a fluidized state, the latter is shown to be virtu- 
ally the only relevant mechanism, and the pressure drop 
follows a simple sinusoidal curve during a cycle (Figs, d 
and |HJ , which corresponds to the weight of the bed per 
unit cross sectional area in the vibrated frame. As a 
consequence, the pressure drop averaged over a cycle in 
the fluidized state is simply the offset of this sinusoidal 
curve, which equals the weight of the bed per unit cross 
sectional area, as in non-vibrated beds. This relation can 
be readily understood by a simple force balance argument 
[Eq. i|17|l]. In a bubbling bed (which appears as 1D-TW 
in our study) , it is during the transient time interval with 
large enough pressure pulsation when the increased ten- 
sile stress breaks up cohesive assembly into agglomerates 

(Fig. CD. 

Note that the compressibility of the gas phase was ig- 
nored in the present study; it would be interesting to 
investigate resonance effects which may arise in vibrated 
beds by allowing for gas compressibility. It was also as- 
sumed that the gas superficial velocity is independent 
of time throughout the oscillation cycle, and this corre- 
sponds to the limit of very large resistance for gas flow 
through the distributor; it would also be interesting to 
examine the case of finite distributor resistance, where 
the superficial gas velocity can be expected to vary cycli- 
cally with plate vibration. 

While the present study has yielded physical under- 
standing of the pressure fluctuations induced by the plate 
and the tensile stress in the particle assembly, an an- 
alytical relation between vibration parameters and the 
agglomerate size is still elusive. Furthermore, we have 
considered only beds of narrow cross sectional areas, and 
assumed the volume-averaged gas phase to be ID. Con- 
sequently, some generic behavior in real vibrated gas flu- 
idized beds, such as horizontal sloshing motion of the 
particles, meandering gas flows around agglomerates, and 
more complicated bubble dynamics, were not allowed to 
occur. By avoiding such complexity, we were able to 
bring out certain basic physics of the bed dynamics. Fu- 
ture studies should investigate the dynamics of higher 
dimensional beds; however, simulation of fluidized beds 
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Nomenclatures 

Hamaker constant 

cohesive Bond number between particles 

cohesive Bond number between particles and base plate 

drag coefficient 

pressure drop across the bed 

grid spacing for coarse-grained variables 

particle diameter 

normal coefficient of restitution 

vibration frequency 

interaction force due to contact 

cohesive force due to van der Waals force 

effective gravitational acceleration 

halo function 

coordination number 

spring stiffness in the normal direction 

spring stiffness in the tangential direction 

mass of individual particle 

solid phase number density 

atmospheric pressure 

gas phase pressure 

Reynolds number based on particle size 
particle radius 
inter-surface distance 
Stokes number 
granular temperature 
bed weight per unit cross sectional area 
volume-averaged gas phase velocity 
volume-averaged solid phase velocity 
superficial gas flow velocity 
minimum fluidization velocity 
relative velocity in normal direction 
relative velocity in tangential direction 
velocity of individual particle 
volume of individual particle 
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Greek symbols 



p 


= 


inter-phase momentum transfer 


S, 5* 


= 


minimum separation distance for cohesion 


A s 


= 


tangential displacement from initial contact 


A„ 


= 


amount of overlap in the normal direction 


e 


= 


porosity (gas phase volume fraction) 


<t> 


= 


solid phase volume fraction 


In 


= 


damping coefficient for dashpot 


r 


= 


vibration intensity 


M 


= 


coefficient of friction 


Pg 




gas phase viscosity 






Poisson's ratio 


Pa 




gas phase mass density 


Ps 




solid phase mass density 


a 




solid phase stress tensor 






tensile strength of the material 


Opiate 




vertical stress on the plate 



TABLE I: Parameter values used. 





Typical values for 
dimensional quantities 


Nominal values for 
dimensionless parameters 


Pa 


1.8xl0 -4 g/(m-s) 




9 


981 cm/s 2 




dp 


50 /im 




Ps 


0.90 g/cm 3 




\fgd~p~ 


2.2 cm/s 




\/d P /g 


2.3xl0~ 3 s 




S 


0.4 nm (= 5*) 


8.0xl0~ 6 


At 


5.6xl0~ 7 s 


2.5xl0~ 4 


Az/dp 




1.5 


e 




0.9 






0.1 






2.0xl0 5 - 2.0xl0 6 


Vp 




0.3 


r 




- 10 


/ 


<~ 100 Hz 


- 0.25 
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